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' ■ ABSTRACT 

Context. Primordial molecules were formed during the Dark Ages, i.e. the time between recombination and reionization in the early 
Universe. They were the constitutents of the first proto-stellar clouds. Standard Big Bang nucleosynthesis predicts the abundances of 
hydrogen, helium, lithium, beryllium, and their isotopes in the early Universe. Heavier nuclei such as carbon, nitrogen, or oxygen are 
formed in trace amounts only. In non-standard Big Bang nucleosynthesis models, it is possible to synthesize larger quantities of these 
heavier elements. The latter are interesting because they can form molecules with a high electric dipole moment which can increase 
the cooling in collapsing protostellar structures. 
(~ l , Aims. The purpose of this article is to analyze the formation of primordial molecules based on heavy elements during the Dark Ages, 

^ I' with elemental abundances taken from different nucleosynthesis models. 

I ' Methods. We present calculations of the full non-linear equation set governing the primordial chemistry. We consider the evolution 

O ' of 45 chemical species and use an implicit multistep method of variable order of precision with an adaptive stepsize control. 

Results. For the first time the cosmological recombination of heavy elements is presented. We find that the most abundant Dark 
^ . Ages molecules based on heavy elements are CH and OH. When considering initial conditions given by the standard Big Bang 

^ I ' nucleosynthesis model, we obtain relative abundances [CH] = «cH/«b = 6.2 x 10"^' and [OH] = ;joH/«b = 1-2 x 10"-' at z = 10, 

where n\, is the total number density. But non-standard nucleosynthesis can lead to higher heavy element abundances while still 
satisfying the observed primordial light abundances. In that case, we show that the abundances of molecular species based on C, N, 
O and F can be enhanced by two orders of magnitude, leading to a CH relative abundance higher than that of HD^ or H2D^. 

^ _ ^ . Key words. Astrochemistry - Early Universe - Cosmology: theory 
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O , /-I'll • 111 1-1 in^, v^iitu^^^u L/oiLn^it; it^av^Liwiia ait iiit^iin^it^iiL uut^ nj lih^ iia- 

The formation of light elements is a central problem of modern ■ Coulomb baiiler and the small density. Therefore, only tiny 

a\ cosmology. Big Bang nucleosynthesis (BBN) provides an im- ^^^^^^^ ^^^^.^^ ^.^ juced in standard BBN, if any. 

portant testing ground for the physics of the primordial Universe. However, a large number of extensions to the standard BBN 

> . The proceedings of nucleosynthesis are relatively well under- j^^^^ ^^^^ suggested, among which we find Non-Standard Big 

stood. Nucleosynthesis is connected with two important cosmo- ^ Nucleosynthesis (hereafter NSBBN) scenarios consider- 

logical events. The first one is the freeze-out, at T ^ 0.8 Me V, ■ ^^^.^^^i^^, j^e early Universe. Some of these 

^ of the weak interactions which interconvert neutrons and pro- ^^^^j^ ^^^^ increased abu ndances of heavy elements (see e.g . 

tons and thus setting the initial proton-to-neutron ratio for nu- j^^j^^^ ^ ^^j^, ^993. jedamz ik et alJl99ltlRauscher et al.[T99l 

cleosynthesis. The second event is the annihilation of thermal ^^j^^^^ ^ ^^^^^.^ 2008 and references therein). The original 

e -e+ pairs m the temperature range 1 MeV - 20 keV, delay- f u„...j; „„.„„u.„; u„..,„„ j„„ 



1. Introduction cleus with the largest binding energy, which is "^He. Beyond 

"^He, charged particle reactions are inefficient due to the ris- 



>< 



pairs m the temperature range 1 MeV - 20 keV delay- j^^tivation for such studies was to obtain an average baryon den 

mg the onset of nucleosynthesis by releasing additional heat but ^. j ^j^^ ^^.^.^^^ ^ ^^^^^^^ 

eliminating the possibility of positron captures. multaneously satisfy the constraints given by primordial abun- 

Thereafter, thermal fusion reactions in the plasma first pro- ^^^^^ observations. While this was shown to be unfeasible (and 
fiiicG fiGiitpniim tiiicIgi* 

unnecessary), density fluctuations are present in the cosmic mi- 

p + ji— >D + y (1) crowave background. The density fluctuations affecting BBN are 

of different magnitude, though, but many different mechanisms 

and then helium (^He, "^He), lithium (^Li) and beryllium (^Be). have been suggested, acting during the inflationary phase, during 

Cosmological expansion plays a crucial role during the whole baryogenesis or later. Thus, NSBBN with density fluctuations 

process and determines the timescale for the nucleosynthesis remains an interesting model for obtaining primordial heavy el- 

during the cooling of the Universe. Observations agree well with ements. 

the pred icted abundances of the standard BBN model. See for During the recombination period the nuclei became progres- 

example iBurles et al.l ( 1200 ll) . but see also ISteigmanI (l2007h for sively neutralized, which led to molecular formation. However, 

the latest observations (and in particular his Fig. 13). at early epochs, where a total absence of dust grains appears 

Due to the proton-rich freeze-out from high temperature, justified, the chemistry is different from the typical interstellar 

standard BBN predicts the highest abundance for the nu- medium astrochemistry. 
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Several groups proposed an a ssembled, comprehen sive set of 
reactions for the early Universe (Lepp & Shull 1984; 'Puv et al.' 
1993; Gain & Palla 1998; Standi et al. 1998b; Signore & Puv 
1999; Lepp et al. 2002; Pfenniger & Puv 2003; Pu v & Signord 
2007; .Glover & AbeL2008; .Glover & Savin, 2009.) . The chemi- 
cal network is coupled with the matter and radiation tempera- 
tures and with the matter density. Primordial chemistry of heavy 
nuclei has been poorly studied, although Lipovkaet al. (2007) 
analyzed the possibility to d etect rotational lines of primor- 
dial CH, and iPuv et al] (|2007|) analyzed the possibility to form 
molecular fluorine HF in the early Universe. 

In this paper we will discuss the chemistry of eight elements 
(hydrogen, deuterium, helium, lithium, carbon, nitrogen, oxy- 
gen and fluorine) in the early Universe, for two different con- 
texts: the standard Big Bang nucleosynthesis model and non- 
standard nucleosynthesis based on baryon density fluctuations. 
Heavier elements are particularly interesting because they can 
form molecules with a high electric dipole moment. Even if their 
abundances are proven to be quite low, they could still furnish an 
interesting cooling agent in the gravitational collapse of proto- 
structures, as molecular abundances are prone to drastic varia- 
tions in such environments. In Sect. |2] we summarize the stan- 
dard and non-standard BBN models we used and the predicted 
abundances. We describe and calculate, in Sect. [3] the mecha- 
nisms of recombination and primordial chemistry which act for 
the eight elements H, D, He, Li, C, N, O and F. In Sect. |4] we 
discuss our results. 



where is the number density of species f and Mb the total num- 
ber density. These values are the initial conditions for our calcu- 
lations of the Standard Big Bang Chemistry (hereafter SBBC). 
The relative abundances for the light nuclei H, D a nd He are in 
accordance with the observational constr aints (see Buries etaTI 
120011; ICvburt et al.ll2003l; ISteigmanllloOTi) . while lithium seems 
to be more abundant than deduced from observations. Indeed, 
Steigman (2007) adopts the following primordial abundance^ 
(D/H)p = 2.68+^■gxl0-^ Yp = 0.240 + 0.006 (< 0.251 +0.002) 
and 12 + log(Li/H)p = 2.1 ± 0.1 (2.5 + 0.1), while the values 
of first column of Table [1] give, after a straightforward calcula- 
tion: (D/H)p = 2.35 X 10-\ Yp = 0.246 and 12 + log(Li/H)p = 
3.3. Let us note that lithium has been observed only in the ab- 
sorption spectra of very old, very metal-poor stars (Population 
II stars); they are considered as ideal for probing the primordial 
abundance of lithium. But lithium is a fragile nucleus which is 
easily destroyed in the interiors of stars. Therefore, if one ad- 
mits a large primordial value for lithium, one must also assume 
a large and uniforme depletion in stars, over a range of stellar 
masses. 

In Table [T] we also quote the amounts of the elements C to F 
produced in this standard BBN, as they are included in the SBBC 
calculations. 

Table 1. Relative abundances [^] - n^/rih of the elements at the 
end of Big Bang nucleosynthesis for the SBBN and for two non- 
standard BBN scenarios. 



2. Big Bang nucleosynthesis 

2.1. Standard model of primordial nucleosynthesis 

Standard BBN assumes a homogeneous baryon density through- 
out the expanding and cooling Universe. The abundance ratio of 
protons and neutrons is set by the freeze-out of the weak in- 
teraction and subsequent decay of the neutrons until the onset 
of nucleosynthesis. At initially high temperature a nuclear sta- 
tistical equilibrium is established, favoring the formation of the 
strongly bound "^He as soon as the photodisintegrations cease. 
All neutrons are consumed but a large fraction of protons re- 
mains because of the higher initial proton abundance. Due to the 
high entropy, i.e. low density, the triple-or reaction is inefficient 
and the formation of elements beyond He is suppressed. For this 
reason, there are few accounts of standard BBN heavy element 
abundances in literature. 

Here, we will consider BBN results bv lWinteled(l2007h . ob- 
tained with a code based on a modified version of the Basel 
network code. Modifications include a consistent evolution of 
temperature and density in the ea rly unive r se, in cluding weak 
freeze-out, based on methods by iKawanol ( |1992|) . and an im- 
proved r eaction network with updated reaction rates, for de- 
tails see IWinteleJ (12007 '). The abundances were calculated for 
77 - 6.22 X 10"'", in accordance with the latest result of the 
Wilkinson Microwave Anisotropy Probe (WMAP) experiment 
(flb/i^ = 0.02273 ± 0.00062, giving 7710 = 6.225 + 0.170, see 
iDunklev et al.ll200"9l) . The weak freeze-out process was carefully 
followed until the nucleosynthesis phase, thus setting the initial 
neutron and proton abundances consistently. The first column of 
Table[T]gives the relative abundances for this SBBN case. Here 
and throughout this paper, we consider relative abundances for 
species ^ such as: 



SBBN NSBBN 







(/v = 0.8; R = 10) 


(A 


= 10-^ R = 


1000) 


[H] 


0.889 


0.889 




0.888 




[D] 


2.092 X 10-^ 


2.45 X 10-5 




2.13 X 10 


-5 


[He] 


0.111 


0.111 




0.112 




[Li] 


1.77 X 10-' 


1.97 X 10-' 




2.30 X 10 


-9 


[C] 


2.51 X lO^''' 


4.00 X 10-'^ 




8.45 X 10- 


14 


[N] 


2.32 X 10^'* 


2.46 X 10-"^ 




3.44 X 10- 


14 


[O] 


3.22 X 10-" 


3.37 X 10-" 




8.20 X 10- 


17 


[F] 


3.28 X 10"" 


3.61 X 10--'' 




1.63 X 10- 


24 



2.2. Non-Standard Big Bang Nucleosynthesis 

While heavy elements are not produced at high levels in 
SBBN, BBN models assuming density fluctuations allow for 
different baryon densities pi,2,... in different zones, leading to 
altered nucleosynthesis. A large number of possibilities for 
creating small scale density perturbations in the very early 
Univer se have be en suggested in literat ure (lAffleck & Dind 
198 5 |ABBlegate'& Hogan 1985; Mala nev & Fowled Il988l:_ 
Dolgov & Silkll993i; iMatsuura et al.i2004trKhlo pov & KouvarisI 
2008). Here, we are not focussing on a specific origin but just 
assume the occurrence of such fluctuations and use the geom- 
etry as open parameter. As customary, we apply a two-zone 
model where the densities pi and p2 of the zones are given by 
the density ratio R - p\lp2 - rf^^ lTf^\ the volume fraction 
< /v < 1 of zone 1, and the additional constraint that the 
averaged density has to reproduce the WMAP value (as in the 
SBBN): 77*'^^P = 6.22 = //("/v + '/'^'(l - /v)- This leaves two 
open parameters, R and /v, but we are further limited by the ob- 
served primordial light abundances. Assuming that the Li abun- 
dance is only a weak constraint because of the complicated stel- 



(2) 



Alternative abundances for "^He and 'Li are written in parentheses. 
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lar depletion of Li, several regions in the /?,/v-space remain al- 
lowed. 

To study the impact of elevat e d leve ls of heavy nuclei, we 
used the same code by IWinteleil (|2007|) as for the SBBN but 
followed nucleosynthesis in zones promising for synthesizing 
heavy elements. Table [1] shows the results for two models, rep- 
resentative for two extreme cases. The abundances shown in the 
table are akeady the final abundances, mixed from the two zones 
in each model: 



/v/;W[^](i)+(l_/^)^(2 ) 

„WMAP 
''10 



(2) 



(3) 



The fi rst case (fy - . 8, R = 10) is similar to the scenario 
studied in iRauscher et alj (Il994l) . The density of pi has to stay 
close to the global tj^J ^^, while p2 is ten times lower. Following 
IRauscher et al.l (1 19941) . to maximize the production of heavy nu- 
clei we assume complete diffusion of the uncharged neutrons out 
of the high-density region. The electrically charged protons re- 
main trapped in that region because of their much shorter mean 
free path. Therefore, baryon density inhomogeneities become lo- 
cal variations of the neutron-to-proton ratio, with free neutrons 
left over in the low-density zone after the initial formation of 
"^He. Neutron captures can then produce neutron-rich isotopes, 
allowing to bypass the slow triple-a reaction and leading to 
heavy elements. Our results are consistent with previous results 
from literature. When obtaining H, D, "^He abundances close to 
the observed ones, the ^Li abundance becomes even higher than 
in the standard model. At the same time, the heavy element abun- 
dances remain at low levels, only very slightly higher than for the 
SBBN. This was already pointed out in Rau scher et al. (1994), 
where it was also found that the light element constraints prevent 
considerable formation of heavy nuclei. 

The second case investigated here (fy - 10"^, R - 1000) 
comprises the other end of the allowed spectrum, tiny pockets 
of extremely high density embedded in a background with al- 
most standard density. This is simila r to the scenario originally 



introduc ed bv | Jeda mzik et aL 
[MatsuuraiLalJ (|2004, 20051 



I and later also discussed by 
In the high density pocket 



the path to heavy elements is opened by an efficient triple-o' re- 
action. At the same time, increased destruction of ^Li keeps its 
abundance low. This simultaneously allows to achieve light el- 
ement abundances compatible with those of the SBBN and to 
increase the heavy element production considerably. However, 
due to the small volume fraction of the high density zone (which 
is necessary to obtain the proper global baryon density) the final 
heavy element abundances after complete mixing remain small, 
although two orders of magnitude higher than for the standard 
case. 

In the following, we will use the abundance values from both 
NSBBN scenarios as initial conditions for our calculations of 
the Non-Standard Big Bang Chemistry (hereafter NSBBCl and 
NSBBC2 respectively). 



3. Dark Ages chemistry 

3.1. Equations of evolution 

The primordial gas is a mixture of hydrogen, deuterium, helium, 
lithium, ... so there are many possibilities of reactions. Change 
in the number density of a given chemical species due to chem- 
ical reactions depends on the densities of the species involved 
in these reactions and on the reaction rates, which themselves 
depend on the matter and radiation temperatures. But of course 



cosmic expansion also plays a crucial role. Thus, it is necessary 
to take into account the following set of differential equations 
(see for example Puv et al. 1993) in the context of the expanding 
Universe, characterized by the scale factor a: 



dt 

dt 

d«b 
dt 



1 da 

— :r^^ 

a dt 

^ J. ^ m + T ^ r W r ^ m) ^ m 

a at 3 nif^c n\, n\, at 

1 dfl 



dnb _ Ida ^ IdnA 



dn^ _ ^ 1 dfl ^ / dn^ 
dt ~ a dt"^ \ dt 



(4) 
(5) 

(6) 
(7) 



chem 



In these equations, is the radiation temperature, the gas 
temperature, n\, the total number density and the number den- 
sity of species ^. The right side of equation (|4| and the first 
term of the right side of equation (|5]l represent the decrease of 
Ti and due to the expansion. The second term of the right 
hand side of ^ is the energy transfer from radiation to the 
gas vi a Compton diffiision of th e CMB photons on the elec- 
trons dKompaneetsll 19571 iPeeblesi, 1968) . In this term, crj is the 
Thomson cross section, the radiation constant, We the elec- 
tronic mass, c the speed of light and «e the electron number den- 
sity. Note that we do not consider the energy transfer between 
gas and radiation via molecular heating and cooling functions, 
since it has been shown that this contribution to the evolution of 
the gas temperature is negligible ( iPuv & Signorell996lll997h . In 
equations ^ and (Q, the first term of the right hand side char- 
acterizes again the density decrease due to cosmic expansion, 
while the second term is the contribution of chemical reactions: 



dn^ 
"dT 



chem 



(8) 



A;^,^, is the rate of the reaction with reactants and ^2- 

This system of ordinary differential equations governing 
the chemical abundances is stiff. We use an implicit multi- 
step method of variable order of precision with an adaptive 
stepsize control. This method has excellent stability proper- 
ties and is widely used for solving chemical kinetic problems 
(Hindmarsh & Petzold 1995). We solve the set of equation (HJi 
to (|7]i from the initial redshift Zinit - 10'^, when the Universe 
was still totally ionized, to the final redshift Zend = 10 (about the 
epoch of reionization). 

Moreover, we consider the standard ACDM model, with 
Hubble parameter //q - 71 km s"' Mpc"', total matter den- 
sity Qm = 0.27 (including dark matter density ^dm - 0.226 and 
baryon density fib = 0.044) and dark energy density Qa = 0.73 
(iKomatsu et al.ll2009l) . 

3.2. Cosmological recombination 

Primordial chemistry begins with the appearance of the first neu- 
tral species. Once a neutral species is formed, charge transfers 
with ions become efficient and lead to the formation of other 
neutral species and then molecular ions. Figure[T]shows the evo- 
lution of the chemical abundances during the successive peri- 
ods of recombination for the light elements H, D, He and Li 
in the standard model. Helium nuclei are the first to recombine 
because of their high electronic binding energjQ. The recom- 



- New re sults on the compl ex helium and hydrogen recombinations 
are given in lFendt et alj ( l2009h (see also references therein). 



4 



P. Vonlanthen et al.: Chemistry of heavy elements in the Dark Ages 



He He 



10000 



1000 10000 
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Fig. 1. Successive cosmological recombination for hydrogen, 
deuterium, helium and lithium in the standard model of primor- 
dial chemistry. 



bination and photoionization rates for H and D are taken from 
[Abel et al. ( 1997). We use the rates of GalH & Palla ( 1998) for 
He and for the reactions Li^ + e Li + y. Higher lithium re- 
combinations and photoionisations are treated with the recombi- 
nation rate s oflVern er et al. ( 1996) and the photoionisation cross 
sections of V erner & Ferlandi (11996). The recombination pattern 
for the two non-standard cases NSBBCl and NSBBC2 are very 
similar, since the initial abundances of H, D, He and Li are al- 
most identical to the SBBC model. We define the recombina- 
tion redshift Zrec as the redshift when the abundance of a given 
species equals the one of its corresponding ion. Table |2] summa- 
rizes the different recombination redshifts for the standard cal- 
culation. Due to its low ionization energy and to the fact that 
charge transfer reactions between lithium and hydrogen remain 
active, lithium does not recombine totally, and both neutral and 
ionized Li tend to match their abundances (see Fig.[T]i. 



Table 2. 

helium. 



Recombination redshifts for hydrogen, deuterium and 



SBBC 


Zrec 


(He+) 


6102 


Zrec 


(He) 


2604 


Zrec 


(H) = Z,ee (D) 


1425 



Moreover, we also investigated, for the first time, the pri- 
mordial recombinations of the carbon, nitrogen and oxygen ions. 
The photoion ization rates are calculated from the cross sections 
evaluated bv iVerner et al.l fl996). The recombination rates are 
taken from|\WL& Ferland ( 1996) for the CVI, CV, CIV, NVII, 
NVI, NV, OVIII OVII and OVI recombinations. They come 
from the work of iPequignot etaTI ( fT99l for the CIH, CII, NIV, 
NIII, Nil and for all other oxygen recombinations. Finally, the 
recombination rates come from ,Omukai (2000) for the CI re- 
combination and from the umist databasj] ( Woodall et al.ll2007l) 
for the NI recombination. We plot in Figs.|2][3]and|4]the succes- 
sive recombinations of carbon, nitrogen and oxygen ions in the 
SBBC model. Recombinations in the NSBBC scenarios are very 
close to the standard recombinations, the only difference being 



Fig. 2. Successive recombinations of the carbon ions in the stan- 
dard model. 
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Fig. 3. Successive recombinations of the nitrogen ions in the 
standard model. 
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Fig. 4. Successive recombinations of the oxygen ions in the stan- 
dard model. 



a shift in relative abundances corresponding to the difference in 
the initial conditions. 



p. Vonlanthen et al.: Chemistry of heavy elements in the Dark Ages 



5 



Since the ionization potential of carbon (Ic - 1 1 -26 eV) is 
smaller than that of hydrogen (/h=13.6 eV), recombination 
occurs later than recombination. Figure |2] shows that carbon 
is singly ionized at the beginning of the primordial chemistry. 
We conclude that it is not necessary to compute all the succes- 
sive recombinations when calculating the evolution of the car- 
bon species. Thus, we will consider that all the carbon ions are 
present as at zinit - 10'*. For similar arguments drawn from 
Figs. |3] and IH we will start our calculations of the primordial 
chemistry of nitrogen and oxygen with the singly ionized nitro- 
gen and oxygen ions at Zinit. 



3.3. Molecular formation 

After hydrogen recombination, the decrease of the electron 
density leads to the inefficiency of thermal coupling between 
matter and radiation. The chemical processes occuring during 
this post-recombination period are essentially coUisional (ion- 
izations, radiative recombinations, attachments, etc.) and ra- 
diative, due to the active presence of the CMB photons. The 
chemical network leads gradually to the formation of the first 
molecules in the Universe (Lepp & Shull 1984; Puv et al. 1993; 
Stancil et al.li 1998b; Galli & Palla 1998; Signore & Puv 1999; 
Lepp et al.l 12002; Pfenniger & Puv ,2003; Puv & Signore..2007i: 



Glover & Abeill2008HGlover & Savinll2009ir 

In this section, we will consider the primordial chemistry of 
hydrogen, deuterium, helium, lithium, carbon, nitrogen, oxygen 
and fluorine in the SBBN and NSBBN contexts. We list in the 
Appendix the sets of reactions that we consider. 

3.3.1. Hydrogen 

The cosmological production of H2 molecules proceeds mainly 
through two mechanisrns which are catalyzed by and Hj 
dSaslaw & Zipovlll967l; [Peebles & Dickelll968l) . It is generally 
considered that the Hj mechanism produces molecular hydro- 
gen at z ~ 300, while the mechanism is responsible for the 
H2 peak at z ~ 100. 

The react i on rate s we consider for hydrogen are taken from 
iGalli & Pallal (Il998h . except for the rec ombination and pho- 
toionization rates HI and H2 (lAbel et al.iri997i) . and for the re- 
action H4 (H" -H 7 — > H -H e~). F or that reaction, we consider the 
important correction described in lHirata & Padmana bhan ( 2006) 
that takes into account the effects of the non-thermal radiation 
emitted during hydrogen recombination. 

Figure |5] shows the evolution of the hydrogen chemistry in 
the standard model. The amount of H2 molecules created during 
the Dark Ages is important and converge to the following val- 
ues at z = 10: [H2]sBBC = «H,/w b - 2 . 7 xlO ~'^. The effect of the 
correction of Hirata & Padmanab hanI (I2006h is shown on Fig.|5] 
The enhancement of the rate H4 is responsible for the efficient 
destruction of via photodetachment. As a consequence, the 
contribution of the channel to the final amount of molecular 
hydrogen is strongly reduced at z ~ 100 and the value of the 
standard final abundance is lower than the values usually seen in 
the literature (about a factor 4). Table[3]gives the abundances at 
redshifts z = 1000, z = 100 and z = 10. The initial relative abun- 
dance of hydrogen is very close in our three chemistry models 
SBBC, NSBBCl and NSBBC2 (see Table as is the case for 
deuterium, helium and to some extent lithium. For that reason 
we only show the results of the standard chemistry calculation 
for these four light elements. 




"10000 



Fig. 5. Standard primordial chemistry of hydrogen species. The 
dotted lines for H2, and show the abundances we would 
obtain without the correction o f iHirata & PadmanabhanI (l2006h 
relative to the radiation emitted during the hydrogen recombina- 
tion. 



Table 3. Relative abundances n^/nb of hydrogen species at red- 
shifts z = 1000, z = 100 and z = 10 in the SBBC model. 



SBBC 


Species 


z = 1000 


z= 100 




Z= 10 


[H] 


0.888 


0.889 




0.889 


[H+] 


8.5 X 10-* 


7.6 X 10- 


5 


5.8 X 10-5 


[H] 


5.7 X 10"^' 


2.5 X 10^ 


15 


7.2 X lO-"* 


[H2] 


1.2 X lO-'"" 


1.7 X 10- 


-7 


2.7 X 10-' 


[H2+] 


7.7 xlO"^" 


2.7 X 10" 


-14 


6.7 X lO-"* 


[H3^] 


3.0x10-22 


8.7 X 10" 


17 


8.4 X 10-"* 



3.3.2. Deuterium 

We take the reacti on rates for deuterium chemistry from 
iGalli & Pallal ( Il998l) . except for the rec ombination and pho- 
toionization rates Dl and D2 dAbel et al.iri997 ). which are the 
same as HI and H2 (see Appendix A). Figure|6|and Table|4|show 
the evolution of the deuterium chemistry. The fact that the abun- 
dance of molecular hydrogen is reduced compared to the usual 
value of the literature has a direct implication on the HD quan- 
tity. Indeed, reactions D6 and D8 are the most important for- 
mation channels of HD. These reactions are collisions between 
neutral or charged D with H2. As a consequence, the standard 
final amount of HD molecules is also lower (by a factor of 2) 
than the canonical value and we get [HD]sbbc = 1.3 x 10-'". 



Table 4. Relative abundances n^/n^ of deuterium species at red- 
shifts z = 1000, z = 100 and z = 10 in the SBBC model. 



SBBC 


Species 


z = 1000 


z= 100 




z= 10 




[D] 


2.1 X 10-' 


2.1 X 10^ 


-5 


2.1 X 10- 


-5 


[D+] 


2.0 X 10-** 


1.3 X 10- 


-9 


2.0 X 10- 


-21 


[HD] 


7.4 X 10-"* 


7.7 X 10- 


-11 


1.3 X 10" 


-10 


[HD+] 


1.8 X 10-2* 


6.9 X 10- 


-19 


4.6 X 10- 


-20 


[HoD+] 


2.1 X 10-2' 


3.8 X 10- 


-20 


4.7 X 10- 


-20 
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Fig. 6. Standard primordial chemistry of deuterium species. As 
in Fig.|5] the dotted lines for HD and H2D^ show the abu ndances 
without the correction of iHirata & PadmanabhanI (l2006h relative 
to the radiation emitted during the hydrogen recombination. 





Fig. 8. Standard primordial chemistry of lithium species. 



3.3.4. Lithium 

We take the reaction rates for lithium chemistry from 
Galli &Palla (1998). Figure Hand Table |6] show the evolution 
of the lithium species. As said before, lithium never recombines 
completely, and this is the reason why LiH is less abundant than 
LiH+ at low redshift: [LiH]sBBC = 13 x 10"^^ and [LiH+]sBBC = 
6.7 X 10-'*^ atz= 10. 



Table 6. Relative abundances n^/nb of lithium species at red- 
shifts z = 1000, z = 100 and z = 10 in the SBBC model. 



Fig. 7. Standard primordial chemistry of helium species. 



SBBC 


Species 


z = 1000 




z = 100 




z = 10 




[Li] 


1.4 X 10" 


20 


6.6 X 10 


-10 


8.8 X 10- 


U) 


[Li+] 


1.8 X 10- 


9 


1.1 X 10 


-9 


8.9 X 10- 


10 


[Li-] 


4.4 X 10- 


41 


2.7 X 10- 


-22 


6.6 X 10- 


21 


[LiH] 


9.5 X 10- 


49 


5.4 X 10^ 


-21 


1.3 X 10- 


19 


[LiH+] 


1.5 X 10- 


38 


1.4 X 10- 


-28 


6.7 X 10- 


18 



3.3.3. Helium 

The reaction rates we consider for helium are taken from 
iGalli & Pailal (Il998 l). The results of the standard helium chem- 
istry are shown in Fig. |7] In the first non-standard BBN model 
we consider (Sect. 12.2b . the formation of helium nuclei dur- 
ing the primordial nucleosynthesis is extremely efficient in the 
neutron-rich, low-density regions. But primordial chemistry be- 
gins after the zones with different densities are mixed by baryon 
diffusion and the non-standard initial relative abundance of he- 
lium is identical to the standard one (see Table [T]i. The main 
feature of helium chemistry is the production of molecule HeH^ 
([HeH+lsBBC = 4.6 x lO-'"* at z = 10). Table |5] gives the abun- 
dances at redshifts z - 1000, z = 100 and z = 10. 

Table 5. Relative abundances n^/wb of helium species at red- 
shifts z = 1000, z = 100 and z = 10 in the SBBC model 



SBBC 


Species 


z = 1000 


z = 100 


z= 10 


[He] 


0.111 


0.111 


0.111 


[He+] 


3.0 X 10-2" 


3.3 X 10-25 


9.1 X 10-26 


[HeH+] 


1.4 X 10-22 


7.2 X 10-'^ 


4.6 X 10-'"* 



3.3.5. Carbon 

The presence of neutral or ionized carbon at high redshift is 
important to determinate the cosmic background temperature at 
z ?^ in orde r to compare the result with models. ISongaila et al.l 
(ll994i[T995h used this technique to prove the existence of carbon 
in some high redshift diffuse gas. In the second paper, the detec- 
tion of absorption due to the first level of neutral carbon fine 
structure, in a z = 1.776 cloud toward quasar Q1331-I-170, al- 
lows them to deduce that the background temperature is 7.4 ±0.8 
K at z = 1.776. This result agrees with the theoretical valu e 
(7.58 K). The same method was used bv iMolaro et all ( l2002h . 
who studied the fine structure levels of in a Ly-o; system at 
z = 3.025 toward Q0347-3819. Their result (Tcmb = 12.1!3 2 K) 
agrees with the standard temperature limits. 

A gas phase che mical mod e l of t he ch emistry of CH and 
CH^ was presented bv lDalgarnol d 1 9761) . then lPrasad & HuntressI 
(11980 ) presented a gas phase chemical model of the chemistry 
in interstellar clouds, including the C, N and O species. We con- 
sider here the evolution of the following carbon species: C, C^, 
C , CH, CH^, CH2 and CH, . The reaction rates are taken from 
the UMisT database ( Woodafl et al. 2007), except for five reac- 
tions: reaction C7 ( Lipovka et al.ii2007i) . the two charge transfers 
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C8 and C9 between C, H and thei r respective ion s dStancil et alj 
\199S sh. the recombination C48 ("Omukail I2OO0I) and the pho- 
toion isation C49, whose cross section is taken from lVerner et alJ 
(11996 ) (see Table |A.5] of the Appendix). The main reaction lead- 
ing to the synthesis of CH is clearly the radiative association 
C15 (C + H ^ CH + y), but we also include the neutral re- 
action CI 2 and the two associative detachments CI 7 and C27. 
Among the numerous reactions destroying CH, the destructive 
collisional reactions with atomic hydrogen (reactions CIO and 
Cll), and the charge exchange with H^ (reaction C16) are the 
most important. The molecular ion CH^ is created by the radia- 
tive association C^ -H H ^ CH^ + y (reaction C23), the col- 
lisional reactions C14 and C18, and the charge exchanges C16 
(with H+) and CI 9 (with H+), the main ones being CI 6 and C23. 
Destruction of CH^ occurs mainly by collisions with H (reaction 
C25) and dissociative recombination (reaction C26). 

In the reaction sets for carbon, nitrogen and oxygen 
chemistries (see Appendix A), we do not include photoioniza- 
tions or photodissociations of primordial molecules (e. g. CH^ 
-I- y — > C^ + H, NH2 + y — > NHj + e"). These processes are es- 
pecially important in the interstellar medium, where strong UV 
photons constitute the background radiation, whereas the Dark 
Ages molecules are embedded in a background radiation that is 
made of much softer microwave photons. In order to be sure 
not to neglect important photoprocesses, we made test calcu- 
lations including photoionizations of the neutral molecules and 
photodissociations of all the molecules based on C, N and O. 
We computed the rates of these reactions using the CMB in- 
tensity and the constant cross section cr - 10"'^ cm^. Indeed, 
this cross section gives a good order of m agnitude estimate for 
the photoionization rates (I van Dishoeck|[T988i) and constitutes 
an upper limit for photodissociation processes. The results of 
the chemistry calculations were not changed when considering 
these photoprocesses. 

Table |7] gives the abundances of the carbon species at red- 
shifts z = 1000, z = 100 and z = 10 for the three models SBBC, 
NSBBCl and NSBBC2. Figure |9] shows the evolution of these 
abundances as a function of redshift in the standard case and in 
the NSBBC2 case. We observe that the main difference between 
the standard and NSBBC2 runs is a global increase of every rel- 
ative abundance by a factor ~ 30 in the NSBBC2 run. Except 
for this shift, the evolution of every carbon species is quite sim- 
ilar in the two scenarios. This is not surprising, since the carbon 
chemistry is mainly determined by reactions with much more 
abundant hydrogen species like H, H^, H , H2 or Hj , and these 
species have the same abundances in all our models. The results 
of the NSBBCl run are quite close to the standard one. We also 
note that the final amount of CH is about 20 times smaller than 
the LiH abundance in the standard model. This ratio is about 
7.5 for HD+ and H2D+. But we note that in the NSBBC2 case, 
molecule CH is as abundant as LiH and more abundant than HD^ 
and H2D+ (about 4.5 times). 

We note that the formation of carbon molecules is com- 
paratively as efficient as the formation of H2 or HD. Indeed, 
the ratio of the most abundant carbon molecule CH to C is 
[CH]/[C] ~ 2.6 X 10"^ at z = 10 in our three models. For 
comparison, we have the following ratios fot the most abundant 
molecule based on each one of the light elements: [H2]/[H] ~ 
3 X 10-^ [HD]/[D] ~ 6 X 10-^ [HeH+]/[He] ~ 4 x 10"^^ and 
[LiH+]/[Li]~8xlO-^atz= 10. 

iLipovka et alj ( l2007h investigated the formation of carbon 
molecules but with a higher initial ratio [C]/[H] - 10""^. They 
found a final CH abundance of 10 which means a similar effi- 
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Fig. 9. Primordial chemistry of carbon species for the SBBC (up- 
per panel) and NSBBC2 models (lower panel). The main differ- 
ence between the two frames is a global increase of every relative 
abundance by a factor ~ 30 in the non-standard model. Except 
for this shift, the evolution of the carbon species is similar in the 
two scenarios (see text). The results of the NSBBCl model are 
quite close to the standard one. 



ciency toward molecular formation, even if their chemical model 
is different. 



3.3.6. Nitrogen 

We take the reaction rates for nitrogen chemi stry from the umist 
database, except for the photoionisation N41 dVerner et alj 19961) 
(see Table |AT6] of the Appendix) and study the evolution of N, 
N^, NH, NH^, NH2 and NHj. The main reactions responsible 
for the creation of NH are collisions of nitrogen atoms with H2 
and H~ (reactions N5 and N14 respectively). We consider many 
processes that can lead to the destruction of NH, the most impor- 
tant of them being collisions with H (reaction N4) and the charge 
exchange H+ + NH NH+ + H (reaction NIO). The latter re- 
action is at the same time one of the main processes responsible 
for the synthesis of the molecular ion NH^, together with the two 
collisional reactions N6 (between N and H;!^) and N7 (between 
N+ and molecular hydrogen). Among the difFerents processes 
destroying NH^, the dissociative recombination N13 appears to 
be the most efficient. 
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Table 7. Relative abundances n^jn^ of carbon species at redshifts z - 1000, z = 100 and z = 10 in the SBBC (left), NSBBCl 
(middle) and NSBBC2 (right) models. 





















NSBBCl 










NSBBC2 






Species 


z = 1000 




z= 100 




z= 10 




Z = 1000 




z= 100 




z= 10 




Z = 1000 




z = 100 




z = 10 




[C] 


7.9 X 10- 


re 


2.4 X 10- 


-15 


2.4 X 10- 


-15 


1.3 X 10- 


15 


3.8 X 10- 


15 


3.8 X 10- 


-15 


2.7 X 10- 


14 


8.0 X 10- 


-14 


8.1 X 10- 


-14 


[C+] 


1.7 X 10- 


15 


1.4 X 10- 


-16 


1.1 X 10 


-16 


2.7 X 10- 


15 


2.3 X 10- 


16 


1.7 X 10- 


-16 


5.7 X 10- 


14 


4.8 X 10- 


-15 


3.7 X 10- 


-15 


[C-] 


7.6 X 10- 


35 


1.1 X 10- 


-24 


6.7 X 10- 


-25 


1.2 X 10- 


34 


1.8 X 10- 


24 


1.1 X 10 


-24 


2.6 X 10- 


33 


3.9 X 10- 


-23 


2.3 X 10- 


-23 


[CH] 


1.2 X 10- 


22 


1.2 X 10- 


-21 


6.2 X 10- 


-21 


2.0 X 10- 


22 


2.0 X 10- 


21 


9.8 X 10 


-21 


4.2 X 10- 


21 


4.2 X 10- 


-20 


2.1 X 10- 


-19 


[CH+] 


3.7 X 10- 


23 


3.4 X 10- 


-24 


3.1 X 10 


-24 


5.8 X 10- 


23 


5.5 X 10- 


24 


5.0 X 10- 


-24 


1.2 X 10- 


21 


1.2 X 10 


-22 


1.1 X 10 


-22 


[CH2] 


7.9 X 10- 


35 


6.9 X 10- 


-29 


1.1 X 10 


-28 


1.3 X 10- 


34 


1.1 X 10 


28 


1.7 X 10- 


-28 


2.8 X 10- 


33 


2.3 X 10- 


-27 


3.7 X 10- 


-27 


[CHJ] 


2.7 X lo- 


34 


1.5 X 10- 


-28 


4.5 X 10- 


-29 


4.3 X 10- 


34 


2.5 X 10- 


28 


7.2 X 10- 


-29 


9.2 X 10- 


33 


5.2 X 10- 


-27 


1.5 X 10 


-27 



Figure \W\ shows the evolution of the nitrogen chemistry 
as a function of redshift in the standard and NSBBC2 cases. 
The molecular ion NH^ dominates the nitrogen chemistry from 
z ~ 600 down to z ~ 80, and is then surpassed by NH, which has 
a peak around z = 40. The final relative abundances at z = 10 are 
[NHlsBBC = 3.4 X 10-25 and [NH]nsbbc2 = 5.1 x lO-^l As for 
carbon, the main difference between the SBBC and NSBBC2 
scenarios is a global increase of every relative abundance (by 
a factor about 150) in the NSBBC2. Apart from this shift, the 
evolution of all the nitrogen species is quite similar in the two 
models. The results of the NSBBCl run are almost identical 
to the standard one, since nitrogen nuclei are produced in the 
same amount in the two corresponding nucleosynthesis scenar- 
ios (Table[T]i- Table[8]gives the abundances at redshifts z = 1000, 
z = 100 and z = 10 for the three models. 

The formation of molecules based on N is less effective than 
the formation of carbon molecules. When considering the most 
abundant species, NH, one has indeed: [NH]/[N] ~ 10"^ for the 
three scenarios. 



3.3.7. Oxygen 

We consider the following species for oxygen chemistry: O, O^, 
O , OH, OH+, H2O and H2O+. Table|AJ]shows the reactions we 
include in our c hemical model. The rate s for reactions 06 1 and 
062 come from lLe Bourlot et al.l (Il993h. The rat e for recombi- 
nation 058 is taken from Peauignot et al. ( 199 lb and the cross 
section for photoionisation 059 from Verner et akl ( 119961) . AU 
other reaction rates come from the umist database. The results 
of oxygen chemistry are shown in Fig.fTTIand Table|9]shows the 
relative abundances of the oxygen species at redshift z = 1000, 
z = 100 and z = 10 for the three models. OH and OH^ are the 
main molecular species based on oxygen at the end of the Dark 
Ages. In order to synthesize OH, we consider, among many other 
reactions, the radiative association 09, collisions between H and 
O" (reaction OlO) and between H" and O (reaction 013), and 
the collisional reaction 015. Note that the first two processes 
are clearly the most effective. Numerous reactions destroy OH 
molecules, but collisions with H (reactions Oil, 012 and 047) 
and H^ (reaction 014) dominate. The molecular ion OH^ is pro- 
duced via the charge exchange 014 (H+ + OH -> OH+ + H), the 
other reactions being less important. OH^ is mainly destroyed by 
dissociative recombination (reaction 022). The collisions with 
hydrogen atoms (OH^ -H H ^ O^ -H H2) are efficient at very 
high redshift only (z > 300). The very strong depletion of O" is 
caused by the photoprocess O + y ^ O -H e" (reaction 05). But 
the UMIST rate for that process assumes an interstellar UV back- 
ground which obviously is inappropriate for our study. We thus 
expect a final O" abundance much higher than plotted in Fig.fTTl 
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Fig. 10. Primordial chemistry of nitrogen species for the SBBC 
(upperpanel) and NSBBC2 models (lower panel). The formation 
of molecules based on nitrogen is less effective than the forma- 
tion of carbon molecules. The main difference between the two 
models shown here is a global increase of every relative abun- 
dance by more than two orders of magnitude in the non-standard 
model. Except for this shift, the evolution of the nitrogen species 
is similar in the two scenarios. The results of the NSBBCl model 
are identical to the standard ones. 



The tendency toward molecular formation is quite strong: we 
observe [0H]/[0] ~ 4 x 10"^ at z = 10 in our three standard and 
non-standard calculations, and even if the initial oxygen abun- 
dances at z = lO'* are very weak, OH is only one order of magni- 
tude less abundant than HD+ and H2D+ in the NSBBC2 run. We 
note the formation of primordial water molecules around z - 30. 
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Table 8. Relative abundances n^/wb of nitrogen species at redshifts z - 1000, z - 100 and z = 10 in the SBBC (left), NSBBCl 
(middle) and NSBBC2 (right) models. 









SBBC 












NSBBCl 










NSBBC2 






Species 


z = 1000 




z= 100 




z= 10 




Z = 1000 




z= 100 




Z = 10 




Z = 1000 




z = 100 




z= 10 




[N] 


2.3 X 10- 


16 


2.3 X 10- 


-16 


2.3 X 10- 


-16 


2.5 X 10- 


16 


2.5 X 10- 


16 


2.5 X 10- 


-16 


3.4 X 10- 


14 


3.4 X 10- 


14 


3.4 X 10- 


14 


[N+] 


1.1 X 10- 


20 


8.1 X 10- 


-22 


6.8 X 10- 


-22 


1.1 X 10- 


20 


8.6 X 10- 


22 


7.2 X 10- 


-22 


1.6 X 10- 


18 


1.2 X 10 


19 


1.0 X 10- 


19 


[NH] 


2.9 X 10- 


30 


3.4 X 10- 


-27 


3.4 X 10- 


-25 


3.1 X 10- 


30 


3.6 X 10- 


27 


3.7 X 10- 


-25 


4.4 X 10- 


28 


5.0 X 10 


25 


5.1 X 10 


23 


[NH+] 


5.4 X 10- 


31 


1.6 X 10- 


-26 


2.1 X 10- 


-27 


5.7 X 10- 


31 


1.7 X 10- 


26 


2.3 X 10- 


-27 


8.1 X 10- 


29 


2.4 X 10- 


24 


3.2 X 10- 


25 


[NH2] 


5.8 X 10- 


44 


3.4 X 10- 


-39 


5.2 X 10- 


-35 


6.2 X 10- 


44 


3.6 X 10- 


39 


5.6 X 10- 


-35 


8.8 X 10- 


42 


5.1 X 10 


37 


7.9 X 10- 


33 


[NH+] 


2.1 X 10- 


42 


1.1 X 10- 


-31 


5.3 X 10- 


-33 


2.3 X 10- 


42 


1.2 X 10- 


31 


5.8 X 10- 


-33 


3.3 X 10- 


40 


1.6 X 10 


29 


8.1 X 10- 


31 



Table 9. Relative abundances n^/n^ of oxygen species at redshifts z = 1000, z = 100 and z = 10 in the SBBC (left), NSBBCl 
(middle) and NSBBC2 (right) models. 









SBBC 












NSBBCl 










NSBBC2 






Species 


z = 1000 




z= 100 




z= 10 




z = 1000 




z= 100 




z= 10 




z = 1000 




z= 100 




z= 10 




[0] 


3.2 X 10- 


19 


3.2 X 10- 


-19 


3.2 X 10- 


-19 


3.4 X 10- 


19 


3.4 X 10- 


19 


3.4 X 10- 


-19 


8.2 X 10- 


17 


8.2 X 10- 


-17 


8.2 X 10- 


-17 


[0+] 


2.7 X 10- 


22 


6.0 X 10- 


-24 


0.0 




2.8 X 10- 


22 


6.2 X 10- 


24 


0.0 




6.9 X 10- 


20 


1.5 X 10- 


-21 


0.0 




[0] 


4.1 X 10- 


28 


6.7 X 10- 


-32 


6.4 X 10- 


-35 


4.2 X 10- 


28 


7.0 X 10- 


32 


7.0 X 10- 


-35 


1.0 X 10- 


25 


1.7 X 10- 


-29 


1.7 X 10- 


-32 


[OH] 


1.8 X 10- 


26 


2.5 X 10- 


-24 


1.2 X 10- 


-23 


1.9 X 10- 


26 


2.6 X 10 


24 


1.3 X 10- 


-23 


4.7 X 10- 


24 


6.2 X 10- 


-22 


3.1 X 10 


-21 


[0H+] 


1.3 X 10 


28 


8.8 X 10- 


-26 


5.2 X 10- 


-26 


1.4 X 10- 


28 


9.2 X 10- 


26 


5.4 X 10- 


-26 


3.3 X 10- 


26 


2.2 X 10- 


-23 


1.3 X 10 


-23 


[H2O] 


1.1 X 10 


37 


1.0 X 10- 


-27 


3.3 X 10- 


-28 


1.1 X 10- 


37 


1.1 X 10 


27 


3.7 X 10- 


-28 


2.8 X 10- 


35 


2.7 X 10- 


-25 


9.2 X 10- 


-26 


[H2O+] 


2.9 X 10 


39 


1.1 X 10- 


-29 


4.4 X 10- 


-31 


3.1 X 10" 


39 


1.2 X 10 


29 


4.9 X 10- 


-31 


7.6 X 10- 


37 


2.8 X 10 


-27 


1.2 X 10 


-28 



They are mainly synthesized by reaction 047 (H + OH — » H2O 
+ y) but at lower redshifts their abundance severly drops due to 
the destructive effect of reaction O60 (H+ + H2O -> H2O+ + H). 
Again, as for carbon and nitrogen, the SBBC and NSBBC2 runs 
differ mainly by a global shift in relative abundance (by more 
than two orders of magnitude). 

3.3.8. Fluorine 

Even if the abundance of primary fluorine is much below the 
abundances of primordial C, N and O at the end of the Big Bang 
nucleosynthesis, the molecule HF could be efficiently formed in 
the early Universe for two reasons: the ionization potential of F 
is greater than that of hydrogen, and the binding energy of the 
hydrogen atom in HF is greater than the binding energy of H2. 
The chemistry of fluorine in the post-recombination epoch has 
been studied in iPuv et alJ ( I2007.) . We use here the same set of 
reactions (see Table |A.8] of the Appendix). 

Figure [12] and Table [10] show the evolution of the fluorine 
chemistry as a function of redshift in the standard and NSBBC2 
cases. We clearly see that the fluorine chemistry is comparatively 
more efficient than the CNO chemistries: [HF]/[F] ~ 10"^ in the 
three models. The enhancement factor for HF in the NSBBC2 
model compared to the standard value is 500. 



4. Summary and Conclusion 

Molecules play an important role in the process of gravita- 
tiona l collapse through the r nechanism of thermal instability 
(see 'Field' '1965; 'SiS '1977; 'Fall & Rees' 'l985t lUehara et alj 
[1996; Puy & Signore 1996; Abel et al. 2000). They provide an 
important cooling mechanism for primordial metal-free gas in 
different primordial contexts su ch as small haloes virializing 
at hi^h redshift (e.g. [B arkana & Loebll2001l; iBromm & LarsonI 
[2OO4; Ciardi & Ferrara 2005) or early structure formation 
(Glover & Jappsen 2007; Maio et al. 2007). 

The H2 molecule has long been recognized as the most 
important cooling agent, despite the absence of dipole mo- 



ment, due to its high relative abundance in comparison with 
other molecules. For this reason H2 molecules are a key 
component of many dynamical situations and particularly 
in the context of a gravita tional collapse. However, several 
studie s iTegmark et al. 1997 :IOmukai & Nishilll998l; lAbel et all 
I2000t iRipamonti et al., ,20021; lYoshida et al.1 120061) consider 
H2 as the only molecular coolant of the primordial gas, al- 
though HP can play an imp ortant role at low temperature. 
Indeed, iPuv & Signord (Il997h showed that, if the primordial 
gas cools below 200 K, HD molecules are the main cooling 
agent, despite their low number abundance. These results 
were co nfirmed in different astrophy sical mediums ([Flowei 
2000; Uehara & Inutsuka 2000; JN^ amura & Umemural 
2000; Flower & Pineau des Forets 200lt i Lipovka et al.' l200jF 
Nagakura & Omukai 2005; Shchekinov & Vasiliev 200g; 
Rip amontii |2007|) . Recently, Prieto et al.l (12008 ) analyzed the 
effects of H2, HD and LiH molecules in the cooling of pri- 
mordial gas. Their simulations clearly showed that the gas, at 
low densities, reaches temperatures about 100 K and that the 
main coolant is H2. But at higher densities (n > 10^ cm"^) HD 
dominates and the gas cools well below 100 K. The effects of 
LiH are negligible in all cases. 



We analyzed in this paper the Dark Ages chemistry of eight 
elements (hydrogen, deuterium, helium, lithium, carbon, nitro- 
gen, oxygen and fluorine), with the aim of determining the 
amount of primordial molecules based on elements heavier than 
lithium. Heavier molecules with high electric dipole moments 
could be important, if created in sufficient quantities, during 
the formation of the first structures in the Universe. We consid- 
ered initial elemental abundances taken from two different nu- 
cleosynthesis contexts: the standard Big Bang nucleosynthesis 
model and non-standard nucleosynthesis based on baryon den- 
sity fluctuations. In both contexts, the baryon-to-photon ratio 77 
(the averaged value for the non-standard model) was in accor- 
dance with the latest result of the WMAP experiment (rf^^^^ = 
6.225+0.170), and the resulting abundances agreed with the ob- 
served primordial abundances of deuterium and helium. 
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Table 10. Relative abundances /«b of fluorine species at redshifts z = 1000, z = 100 and z = 10 in the SBBC (left), NSBBCl 
(middle) and NSBBC2 (right) models. 
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Fig. 11. Primordial chemistry of oxygen species for the SBBC 
(upper panel) and NSBBC2 models (lower panel). OH is clearly 
the most abundant oxygen molecule. Primordial water molecules 
show a peak around z = 30 but are then greatly depleted at 
lower redshifts. The main diff'erence between the two models is a 
global increase of every relative abundance by more than two or- 
ders of magnitude in the non-standard NSBBC2 model. Except 
for this shift, the evolution of the oxygen species is similar in the 
two scenarios. The results of the NSBBCl model are very close 
to the standard ones. 



The most abundant Dark Ages molecules based on elements 
heavier than lithium are CH and OH. Chemistry calculation as- 
suming standard BBN initial abundances yields [CH]sbbc - 
6.2 X 10-21 and [OH]sbbc = 1 .2 x lO^^^ at z = 10. NH molecules 
are less abundant; [NH]sbbc = 3.4 x 10"^^ 

We also made two non-standard chemistry calculations, con- 
sidering two specific cases of inhomogenous BBN. In the first 
one, primordial nucleosynthesis abundances for C, N, O and F 
are only very slightly higher than in the standard case. For that 



Fig. 12. Primordial chemistry of fluorine species for the SBBC 
(upper panel) and NSBBC2 models (lower panel). HF is more 
efficiently created by a factor 500 in the NSBBC2. The results 
of the NSBBCl model are almost identical to the standard ones. 



reason, the chemistry of elements heavier than lithium in that 
context is very similar to the standard chemistry. But in the sec- 
ond non-standard BBN case, the final heavy element abundances 
are enhanced by globally two orders of magnitude. In that case, 
the relative abundances of molecular species based on carbon, 
nitrogen, oxygen or fluorine are modified by essentially the same 
factor, and molecules CH are as abundant as LiH and more abun- 
dant than HD+ and H2D+: [CH]nsbbc2 = 2.1 x lO"''' at z = 10. 
The details of the chemistry (e. g. the most important species, 
the most effective reaction s, . . . ) are weakly affected. 

iGlover & Abell (l2008l) investigated the effects of the chem- 
ical rate coefficient uncertainties in the Dark Ages chemistry 
of light elements. They showed that the large uncertainties in 
the associative detachment and mutual neutralization rates have 
an impact on the thermal evolution of the gas. As regards the 
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heavy molecules discussed in the present work, the main reac- 
tions creating or destroying them are affected by uncertainties 
that can come to a factor 2. Moreover, some of the rates included 
in our model are given for a more restricted temperature range 
than the one of our calculations (from z - lO'^toz = 10 the 
Universe cools from T ~ 30000 K down to a few Kelvins). We 
nonetheless assumed that these rates were valid everywhere in 
that range. Concerning oxygen chemistry, the UMIST rate for 
reaction 05 (O" + y O + e ) is largely not appropriate to the 
Dark Ages chemistry. As one of the two main reactions creating 
OH is H + ^ OH + e , the use of a better rate for 05 could 
increase the OH relative abundance obtained in this work. 

Knowledge of primordial molecular abundances of heavy 
molecules could be important in different contexts, cosmologi- 
cally as well as for the formation of the first structures and stars. 
A cosmo logical influence of primordial molecules was sug- 
gested b V iDubrovichI d 1 9771) . The measure of CMB anisotropics 
gives accurate indications on the fluctuation spectrum which led 
to present great structures. In this way, it is possible that primor- 
dial molecular clouds affect the anisotropy spectrum at small 
scales via resonant scattering of CM B photo ns on primor dial 
molecules dMaoh et alj 1 1994 119961: I Schleiche r'etani2008l for 
a review see |Basu||2007| or|Signore & Puy 2009). The first spec- 



tral line surveys searchi ng for primordial signals have been done 
with the Odin satellite dPersson et al.ll2009i) . Their work may be 
considered as a pilot study for the forthcoming observations with 
the Herschel Space Observatory launched on May 14, 2009. One 
of the most promising molecule has been proven to be HeH^. 
Molecules such as CH, NH or OH have electric dipole moments 
of the same order as HeH^. We have shown in this work that 
their abundances are smaller by several orders of magnitude than 
HeH^ abundance. But the fact that CH could be more abundant 
than molecules such as HD^ and H2D^ could open new interest- 
ing perspectives and detailed calculations have still to be done. 
In addition to molecules, heavy atom s and ions can also be used 
as a tool to probe the early Universe. iBasu et aTl d2004l) consid- 
ered the opacity generated by the scattering of CMB photons on 
heavy atoms like carbon, oxygen, silicon or iron (and on their 
respective ions). They found that the Planck HFI detectors will 
be able to get strong constraints (lO"-' - 10""* solar abundance) 
on the mo st important atoms and ions (CNO) in the interval 
z € [5, 301. lHarwit & SpaansI d2003h also examined the possibil- 
ity to detect the imprint left by the presence of carbon or oxygen 
on the CMB during the Dark Ages. 

Numerical simulations of star formation are a very impor- 
tant and popular problem in modern cosmology. If solved suc- 
cessfully, it can add a lot of information to our knowledge of the 
evolution of all kinds of cosmic structures and the Universe as a 
whole. Star formation is determined by complex chemi cal, dy- 
namic al and thermodynamical processes (see McGreer & BrvanI 
|2008|) . each of which plays an important role during some or all 
parts of the whole process, and the impact of cooling by heav- 
ier molecules on the formation of the first (Population III) stars 
is a new open question. Moreover, it is interesting to note that 
the se first stars were th emselves an important source of metallic- 
ity dCooke et al.ll2009l) , and cooling by heavier molecules could 
be even more important for later star formation. 



Acknowledgements. We acknowledge the PNCG (Programme National 
Cosmologie et Galaxies) and the PNPS (Programme National de Physique 
Stellaire) for their financial assistance, and Daniel Pfenniger, Anton Lipovka 
and Dahbia Talbi for fruitful discussions. This work was in part supported by 
the Swiss National Science Foundation (grant 200020-122287). 



References 

Abel, T., Anninos, P., Zhang, Y., & Norman, M. L. 1997, New Astronoiny, 2, 
181 

Abel, T., Bryan, G. L., & Norman, M. L. 2000, ApJ, 540, 39 
Afileck, I. & Dine, M. 1985, Nuclear Physics B, 249, 361 
Applegate, J. H. & Hogan, C. J. 1985, Phys. Rev. D, 31, 3037 
Barkana, R. & Loeb, A. 2001, Phys. Rep., 349, 125 
Basu, K. 2007, New Astronomy Review, 51, 431 

Basu, K., Hernandez-Monteagudo, C., & Sunyaev, R. A. 2004, A&A, 416, 447 

Bromm, V. & Larson, R. B. 2004, ARA&A, 42, 79 

Buries, S., Nollett, K. M., & Turner, M. S. 2001, ApJ, 552, LI 

Ciardi, B. & Fen'ara, A. 2005, Space Science Reviews, 1 16, 625 

Cooke, J., Cooray, A., Chary, R.-R., et al. 2009, ArXiv e-prints 0902.4602 

Cyburt, R. H., Fields, B. D., & Olive, K. A. 2003, Physics Letters B, 567, 227 

Dalgarno, A. 1976, in Atomic processes and applications. P. G. Burke (ed.), 

North-Holland Publ. Co., Amsterdam, Netherlands, p. 109 - 132, 109-132 
Dolgov, A. & Silk, J. 1993, Phys. Rev. D, 47, 4244 
Dubrovich, V. K. 1977, Soviet Astronomy Letters, 3, 128 
Dunkley, J., Komatsu, E., Nolta, M. R., et al. 2009, ApJS, 180, 306 
Fall, S. M. & Rees, M. J. 1985, ApJ, 298, 18 

Fendt, W. A., Chluba, J., Rubifio-Martfn, J. A., & Wandelt, B. D. 2009, ApJS, 

181, 627 
Field, G. B. 1965, ApJ, 142, 531 
Flower, D. R. 2000, MNRAS, 318, 875 

Flower, D. R. & Pineau des Forets, G. 2001, MNRAS, 323, 672 
Galh, D. & Palla, E 1998, A&A, 335, 403 
Glover, S. C. O. & Abel, T. 2008, MNRAS, 388, 1627 
Glover, S. C. O. & Jappsen, A.-K. 2007, ApJ, 666, 1 
Glover, S. C. O. & Savin, D. W. 2009, MNRAS, 393, 91 1 
Harwit, M. & Spaans, M. 2003, ApJ, 589, 53 

Hindmarsh, A. C. & Petzold, L. R. 1995, Computers in Physics, 9, 34 

Hirata, C. M. & Padmanabhan, N. 2006, MNRAS, 372, 1175 

Jedamzik, K., Fuller, G. M., Mathews, G. J., & Kajino, T. 1994, ApJ, 422, 423 

Kawano, L. 1992, Fermilab-Pub-92/004-A 

Khlopov, M. Y. & Kouvaris, C. 2008, Phys. Rev. D, 77, 065002 

Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2009, ApJS, 180, 330 

Kompaneets, A. S. 1957, Sov Phys. JETR 4, 730 

Le Bourlot, J., Pineau Des Forets, G., Rouefl", E., & Flower, D. R. 1993, A&A, 
267, 233 

Lepp, S. & Shull, J. M. 1984, ApJ, 280, 465 

Lepp, S., Standi, P. C, & Dalgarno, A. 2002, Journal of Physics B Atomic 

Molecular Physics, 35, 57 
Lipovka, A., Nuiiez-Lopez, R., & Avila-Reese, V. 2005, MNRAS, 361, 850 
Lipovka, A., Saucedo, J., & Lipovka, N. M. 2007, ArXiv e-prints 0703547 
Maio, U., Dolag, K., Ciardi, B., & Tomatore, L. 2007, MNRAS, 379, 963 
Malaney, R. A. & Fowler, W. A. 1988, ApJ, 333, 14 

Maoli, R., Ferrucci, V., Melchiorri, E, Signore, M., & Tosti, D. 1996, ApJ, 457, 
1 

Maoli, R., Melchiorri, F, & Tosti, D. 1994, ApJ, 425, 372 
Matsuura, S., Dolgov, A. D., & Nagataki, S. 2004, Progress of Theoretical 
Physics, 112, 971 

Matsuura, S., Fujimoto, S.-I., Hashimoto, M.-A., & Sato, K. 2007, Phys. Rev D, 
75, 068302 

Matsuura, S., Fujimoto, S.-L, Nishimura, S., Hashimoto, M.-A., & Sato, K. 2005, 

Phys. Rev. D, 72, 123505 
McGreer, I. D. & Bryan, G. L. 2008, ApJ, 685, 8 

Molaro, P., Levshakov, S. A., Dessauges-Zavadsky, M., & D'Odorico, S. 2002, 

A&A, 381, 64 
Nagakura, T. & Omukai, K. 2005, MNRAS, 364, 1378 

Nakamura, F. & Umemura, M. 2000, in The First Stars, ed. A. Weiss, T. G. Abel, 

& V. Hill, 263 
Omukai, K. 2000, ApJ, 534, 809 
Omukai, K. & Nishi, R. 1998, ApJ, 508, 141 
Peebles, P J. E. 1968, ApJ, 153, 1 
Peebles, P J. E. & Dicke, R. H. 1968, ApJ, 154, 891 
Pequignot, D., Petitjean, P, & Boisson, C. 1991, A&A, 251, 680 
Persson, C. M., Maoh, R., Encrenaz, P., et al. 2009, A&A, submitted 
Pfenniger, D. & Puy, D. 2003, A&A, 398, 447 
Prasad, S. S. & Huntress, Jr, W. T. 1980, ApJS, 43, 1 
Prieto, J. P., Infante, L., & Jimenez, R. 2008, ArXiv e-prints 
Puy, D., Alecian, G., Le Bourlot, J., Leorat, J., & Pineau Des Forets, G. 1993, 

A&A, 267, 337 

Puy, D., Dubrovich, V., Lipovka, A., Talbi, D., & Vonlanthen, P 2007, A&A, 
476, 685 

Puy, D. & Signore, M. 1996, A&A, 305, 371 

Puy, D. & Signore, M. 1997, New Astronomy, 2, 299 

Puy, D. & Signore, M. 2007, New Astronomy Review, 51, 411 



12 



P. Vonlanthen et al.: Chemistry of heavy elements in the Dark Ages 



Rauscher, T., Applegate, J. H., Cowan, J. J., Thielemann, F.-K., & Wiescher, M. 

1994, ApJ, 429, 499 
Ripamonti, E. 2007, MNRAS, 376, 709 

Ripamonti, E., Haardt, R, Feirara, A., & Colpi, M. 2002, MNRAS, 334, 401 

Saslaw, W. C. & Zipoy, D. 1967, Nature, 216, 976 

Schleicher, D. R. G., Galli, D., Palla, R, et al. 2008, A&A, 490, 521 

Shchekinov, Y. A. & Vasiliev, E. O. 2006, MNRAS, 368, 454 

Signore, M. & Puy, D. 1999, New Astronomy Review, 43, 185 

Signore, M. & Puy, D. 2009, European Physical Journal C, 59, 1 17 

Silk, J. 1977, ApJ, 211,638 

Songaila, A., Cowie, L. L., Vogt, S., et al. 1994, Nature, 371, 43 

Songaila, A., Hu, E. M., & Cowie, L. L. 1995, Nature, 375, 124 

Stancil, P C, Havener, C. C, Krstic, P S., et al. 1998a, ApJ, 502, 1006 

Stancil, P C, Lepp, S., & Dalgarno, A. 1998b, ApJ, 509, 1 

Steigman, G. 2007, Annual Review of Nuclear and Particle Science, 57, 463 

Tegmark, M., Silk, J., Rees, M. J., et al. 1997, ApJ, 474, 1 

Uehara, H. & hiutsuka, S.-i. 2000, ApJ, 531, L91 

Uehara, H., Susa, H., Nishi, R., Yamada, M., & Nakamura, T. 1996, ApJ, 473, 
L95 

van Dishoeck, E. R. 1988, in Rate coefficients in astrochemistry, ed. Dordrecht 

and Norwell, MA, Kluwer Academic Publishers, 49-72 
Vemer, D. A. & Rerland, G. J. 1996, ApJS, 103, 467 

Vemer, D. A., Rerland, G. J., Korista, K. T., & Yakovlev, D. G. 1996, ApJ, 465, 
487 

Winteler, C. 2007, Master thesis. University of Basel, Switzerland 
Woodall, J., Agiindez, M., Mai'kwick-Kemper, A. J., & Millai', T. J. 2007, A&A, 
466, 1197 

Yoshida, N., Omukai, K., Hemquist, L.. & Abel, T. 2006. ApJ, 652, 6 



Appendix A: Reaction rates for the 
post-recombination chemistry 

We list in Tables [ATI to IA.8I the chemical reactions considered 
in our calculations. 



Table A.l. List of reactions for t he hy drogen chemistry. The 
rates are taken from Galli & Pallal (Il998h . exc ept for the recom - 
bination and photoio nization rates HI and H2 ( Abel et aO 19971) . 
and for reaction H4 (iHirata & Padmanabhanl2006E 



Table A.3. Lis t of reactions for the heUum chemistry. The rates 
are taken from lGalh & Pallal (Il998h . 



HI 


H+ 


+ e —>}l + y 


H2 


H + 


7 ^ H+ + 


H3 


H- 


H — > + 7 


H4 


H" 


f 7 — > H + 


H5 




+ H ^ H2 + 


H6 


H" 


f H+ ^ H^ + 


H7 




+ H+ ^ H + H 


H8 


H + H+ ^ H^ + 7 


H9 




+ y^U + U* 


HIO 


H* 


f H ^ H2 + H+ 


Hll 




+ e- ^ H + H 


HI2 


HI, 


f 7 ^ H+ + H+ + e- 


H13 




+ H2 ^ H+ + H 


HIS 


H^ ^ 


h H+ ^ H+ + H 


H16 


H2 


+ ^ H + 


H17 


H2^ 


h -< H + H + 


H18 


H, 


+ y — > HJ + 


HI9 


Ht 


f H ^ Hi + H2 


H20 




+ ^ H, + H 


H21 


H2 ^ 


h H+ ^ H+ + 7 



Table A.2. Li st of reactions for the deuterium chemistry. We take 
the rates from lGalli & Pallal (11998 ^. except for the recombination 
and pho toionization rates Dl and D2, which are the same as HI 
and H2 (lAbel et al.lll997l) . 





Reaction 




Reaction 


DI 


D+ + — > D + 7 


D2 


D + 7 — > D+ + e" 


D3 


D + H+ ^ D+ + H 


D4 


D+ + H ^ D + H+ 


D5 


D + H ^ HD + 7 


D6 


D + H2 H + HD 


D7 


HD+ + H ^ H+ + HD 


D8 


D+ + H2 H+ + HD 


D9 


HD + H ^ H2 + D 


DIO 


HD + H+ ^ H2 + D+ 


DIl 


HD + H+ ^ H2 + H2D+ 


DI2 


D + H+ ^ HD+ + 7 


D13 


D+ + H ^ HD+ + 7 


D14 


HD+ + 7 ^ H + D+ 


D15 


HD+ + 7 -> H+ + D 


D16 


HD+ + e- ^ H + D 


D17 


HD+ + H2 ^ H2D+ + H 


D18 


HD+ + H2 ^ H+ + D 


D19 


D + H+ ^ H2D+ + H 


D20 


H2D+ +e-^ti + H + D 


D21 


H2D+ + e- ^ H2 + D 


D22 


H2D+ + e- ^ HD + H 


D23 


H2D+ + H2 ^ H+ + HD 


D24 


H2D+ + H ^ H+ + D 





Reaction 






Reaction 


Hel 


He"^+ + e" - 


-> He^ + y 


He2 


He* + 7 — > He+* + e" 


He3 


He+ + e" — i 


He + 7 


He4 


He + 7 — > He+ + e" 


He5 


He + H+ ^ 


He+ +H 


He6 


He+ + H ^ He + H+ 


He7 


He + H+ ^ 


HeH+ + 7 


He8 


He + Hi ^ HeH+ + H 


He9 


He+ + H ^ 


HeH+ + 7 


HelO 


HeH+ + H ^ He + HJ 


Hell 


HeH+ + e" 


He + H 


Hel2 


HeH+ + H2 ^ H| + He 


Hel3 


HeH+ + 7 - 


^ He + H+ 


Hel4 


HeH+ + 7 ^ He+ + H 



Table A.4. List of reactions for the lithium chemistry. We con- 
sider the rates discussed in Galli & Palla (1998), except for the 
recombinations Li23 and Li25 (Verneret al. 1 996). and for the 
photoionisations Li24 and Li26 (Verner & Ferl and|[l996l) . 





Reaction 






Reaction 




Lil 


Li* + e" - 


-> Li + 7 


Li2 


Li + 7 — > 


Li+ + e- 


Li3 


Li+ + H- 


Li + H 


Li4 


Li- + H+ 


^Li + H 


Li5 


Li + e" ^ 


Li" + 7 


Li6 


Li" + 7 — 


> Li + e- 


Li7 


Li + H+ ^ Li+ + H 


Li8 


Li + H+ - 


^ Li+ + H + 7 


Li9 


Li + H- - 


^ LiH + e- 


LilO 


Li- + H - 


-> LiH + e- 


Lill 


LiH+ + H 


^ LiH + H+ 
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LiH + H+ 


^ LiH+ + H 


Lil3 


LiH + H - 


^ Li + H2 


Lil4 


Li+ + H - 


^ LiH+ + 7 
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Li + H+ - 


^ LiH+ + 7 


Lil6 


LiH + H+ 


^ LiH+ + H 
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LiH + H+ 


^ Li+ + H2 


Lil8 


LiH+ + e- 


-i Li + H 
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LiH+ + H 


^ Li + Hi 


Li20 


LiH+ + H 


^ Li+ + H2 


Li21 


LiH+ + 7 


^ Li+ + H 


Li22 
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^ Li + H+ 


Li23 


Li-+ + e" 


^ Li+ + 7 


Li24 


Li+ + 7 - 


. Li-* + e- 


Li25 


Li3+ + e- 


^ Li-+ + 7 


Li26 


Li2+ + 7 ^ Li'+ + e- 



Table A.5. List of reactions for the carbon chemistry. The rates 
com e from the umist data base, except for five reactions: reaction 
C7 (iLipovka et al .1120071) . the two charge transfers C8 and C9 
between C, H and their respective ions (Stancil et al. l99Ba), the 
recombination C48 (Omukai 20Q( i) and the photoionis ation C49, 
whose cross section is taken from lVerner et al.l (Il996h . 





Reaction 




Reaction 


Cl 


C+ + H- ^ H + C 


C2 


He+ + C ^ C+ + He 


C3 


H+ + C- ^ C + H 


C4 


He+ + C- ^ C + He 


C5 


c+ + c- ^ C + C 
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C + e- -> C- +7 


C7 


C- + 7 -> C + e- 
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C+ + H ^ C + H+ 


C9 


C + H+ ^ C+ + H 


CIO 


H + CH ^ C + H2 


ClI 


H + CH^C + H + H 
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C + H2 ^ CH + H 
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H2 + CH -t C + H + H2 


CI4 


C+ + H2 ^ CH+ + H 


CIS 


C + H CH + 7 
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H- + C ^ CH + e- 
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He+ + CH ^ CH+ + He 
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CH+ + e- ^ C + H 
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H + C- ^ CH + e- 
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H2 + CH ^ CH2 + H 
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H+ + CH2 ^ CH+ + H2 


C31 


H + CH+ -> CH+ + H2 
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H2 + CH+ -4 CH* + H 
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Hi + CH ^ CHi + H2 
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He+ + CH2 ^ CH+ + He + H 
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H+ + CH2 ^ CHi + H 
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C + CH2 ^ CH + CH 
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Hi + CH2 ^ CHi + H2 
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C+ + CH2 ^ CH+ + C 


C41 


Cifi +e- -4 CH + H 


C42 


CH* +e-^C + H + H 


C43 


CHi + e- ^ C + H2 


C44 


H- + CH ^ CHj + e- 


C45 


H2 + C- ^ CH2 + e- 


C46 


H2 + C ^ CHj + 7 


C47 


H2 + C+ ^ CHi + 7 


C48 


C+ + e- ^ C + 7 


C49 


C + 7 ^ C+ + e- 
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Table A.6. List of reactions for the nitrogen chemistry. We take 
the rates from the umist database, except for the photoionisation 
N41 (Verneretal. 1996^. 





Reaction 




Reaction 


Nl 


H- + N+ ^ N + H 


N2 


C- + N+ ^ N + C 


N3 


CH + N+ ^ N + CH+ 


N4 


H + NH ^ N + H2 


N5 


Ht + N ^ NH + H 


N6 


H| + N ^ NH+ + H 


N7 


H2 + N+ ^ NH+ + H 


N8 


H2 + NH+ ^ N + H+ 


N9 


He+ + NH ^ N+ + He + H 


NIO 


H+ + NH ^ NH+ + H 


Nil 


HJ + NH ^ NH+ + H2 


N12 


N+ + NH ^ NH+ + N 


N13 


NH+ + e" ^ N + H 


N14 


H" + N ^ NH + e" 


N15 


CH + N ^ NH + C 


N16 


N + CH2 ^ NH + CH 


N17 


C + NH+ -> N + CH+ 


N18 


CH + NH+ -> CH+ + N 


N19 


C + NH ^ N + CH 


N20 


H + NH2 ^ NH +"H2 


N21 


H2 + NH ^ NH2 + H 


N22 


NH + NH ^ NH2 + N 


N23 


H+ + NH -» NHi + H 


N24 


H2 + NH+ ^ NH+ + H 


N25 


H| + NH NHt + H2 


N26 


He+ + NH2 ^ N+ + He + H2 


N27 


He+ + NH2 ^ NH+ + He + H 


N28 


NH+ + NH ^ NHJ + N 


N29 


H+ + NH2 ^ NHJ + H 


N30 


HJ + NH2 ^ NHJ + H2 


N31 


N+ + NH2 -4 NH+ + N 


N32 


NH+ +e-^N + H + H 


N33 


NH+ + e" ^ NH + H 


N34 


H" + NH ^ NH2 + e" 


N35 


C + NH2 ^ NH + CH 


N36 


C + NHJ ^ NH + CH+ 


N37 


CH + NHJ ^ NH + CH+ 


N38 


CH + NH+ ^ NH2 + CH+ 


N39 


CH2 + NHJ ^ NH2 + CH+ 


N40 


N* + e" — > N + y 


N41 


N + 7 — > N"^ + e" 







Table A.7. List of reactions for the oxygen chemistry. We 
consider the r ates of the umist database, except for reactions 
061 and 062 dLe Bourlot et al.lil993) . for recombination 058 
jPequignotet al.l Il991h an d for photoi onisation 059, whose 
cross section is taken from (IVerner et al.|[r996) . 





Reaction 




Reaction 


01 


+ e- ^ Q- + y 


02 


H + 0+ ^ + H+ 


03 


H+ + ^ 0+ + H 


04 


H" + 0+ ^ + H 


05 


0" + y -> + e" 


06 


c- + 0+ ^ + C 


07 


CH + 0+ ^ + CH+ 


08 


NH + 0+ ^ + NH+ 


09 


H + ^ OH + y 


010 


H + 0" ^ OH + e- 


Oil 


H + OH^O + H + H 


012 


H + OH ^ + H2 


013 


H + ^ OH + e" 


014 


H+ + OH ^ 0H+ + H 


015 


H2 + -» OH + H 


016 


H2 + OH ^ + H2 + H 


017 


H2 + 0+ ^ 0H+ + H 


018 


Hj + ^ 0H+ + H 


019 


H+ + OH ^ 0H+ + H2 


020 


H+ + ^ 0H+ + H2 


021 


0+ + OH ^ 0H+ + 


022 


0H+ + e" ^ + H 


023 


He+ + OH ^ 0+ + He + H 


024 


C + OH ^ + CH 


025 


C + 0H+ ^ + CH+ 


026 


CH + ^ OH + C 


027 


CH + 0H+ ^ OH + CH+ 


028 


OH + F ^ HF + 


029 


N + OH ^ + NH 


O30 


N+ + OH ^ 0H+ + N 


031 


NH + ^ OH + N 


032 


NH+ + ^ 0H+ + N 


033 


H2O+ +e-^0 + H + H 


034 


H2O+ + e- ^ + H2 


035 


H2O+ + e- ^ OH + H 


036 


+ H2O ^ OH + OH 


037 


0+ + H2O ^ H2O+ + 


038 


H2 + 6- ^ H2O + e- 


039 


H2 + H2O ^ OH + H2 + H 


O40 


H2 + OH ^ H2O + H 


041 


H2 + 0H+ ^ H2O+ + H 


042 


H+ + OH ^ H2O+ + H 


043 


H| + ^ H2O+ + H 


044 


H+ + OH ^ H2O+ + H2 


045 


H + H2O ^ OH + H2 


046 


H V H2O ^ OH + H + H 


047 


H + OH ^ H2O + y 


048 


H" + OH ^ H2O + e" 


049 


He+ + H2O ^ OH + He + H+ 


O50 


He+ + H2O 0H+ + He + H 


051 


OH + OH ^ H2O + 


052 


0H+ + OH H2O+ + 


053 


0H+ + H2O H2O+ + OH 


054 


H2O + F ^ HF + OH 


055 


C + H2O+ ^ OH + CH+ 


056 


NH + OH ^ H2O + N 


057 


NH+ + OH ^ H2O+ + N 


058 


0+ + e- ^ + y 


059 


+ y -> 0+ + e" 


O60 


H+ + H2O ^ H2O+ + H 


061 


0H+ + H 0+ + H2 


062 


H2O+ + H ^ H+ + H2O 


063 


H+ + H2O ^ H2O+ + H2 







Table A.8. List of re actions for the flu orine chemistry. We take 
the rates discussed in lPuy et aP (l2007h . 





Reaction 




Reaction 




Fl 


F + H ^ HF + y 


F2 


FH + y- 


>F + H 


F3 


F + e" -> F" + y 


F4 


F" + y ^ 


F + e- 


F5 


F" + H ^ HF + e" 


F6 


HF + e" 


F" + H 


F7 


F + H2 ^ HF + H 









